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Abstract. The fragmentation functions of the pion with distinction between DZ , -DJ , and are 
studied in the Field-Feynman recursive model, by taking into account the flavor structure in the excitation 
of quark-antiquark pairs by the initial quarks. The obtained analytical results are compatible with available 
empirical results. The framework is also extended to predict the fragmentation functions of the kaon with 



distinction between (z), (z), Df (z), and Df (z). This work gives a significant modification of 
the original model, and the predictions can be tested by future experiments on the fragmentation functions 
of the kaon. 

PACS. 13.87.Fh - 12.38.Lg - 13.60.Le - 14.40.Aq 



1 Introduction 

The productions of mesons are under investigations in 
various processes nowadays, therefore the fragmentation 
functions of mesons are among the useful basic quanti- 
ties in high energy physics studies. A complete knowledge 
of the quark to meson fragmentation functions with de- 
tailed flavor structure is needed and will be more use- 
ful. In previous studies, the fragmentation functions of 
mesons are usually obtained from experimental data on 
e + e~ annihilation process 1 . However, the situation has 
now been changed because the HERMES collaboration 
has published results of charge separated data for 7r ± pro- 
duction on the proton target [2] • By combining the HER- 
MES data on semi-inclusive DIS (SIDIS) ^-production 
with the singlet fragmentation function , which is well 
determined from e + e~ data, Kretzer, Leader, and Chris- 
tova |2| obtained a complete set of fragmentation func- 
tions of the pion: , , and , without specific 
assumption about relation between favored and unfavored 
transitions, such as DJ + (z) = (1 — z)D^ + (z) in g] and 

Df{z) = (1 - z) n Dl + {z) (with n=2,3,4) in 5 . The 
Kretzer-Leader-Christova (KLC) parametrization results 
of the experimental data are: 



Dl (z) = 0.689z 



-1.039 



(1 



\1.241 



Dl (z) = 0.217*- 1,BOl, (l - z) 



\2.037 



(2) 
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Df (z) = 0.164z- 1 - 927 (l - z) 2 ' 886 . (3) 

Although the evolution behaviors can be calculated 
perturbatively, the actual form of the fragmentation func- 
tions is non-perturbative. Therefore the fragmentation func- 
tions arc helpful to elucidate the fundamental features of 
hadronization. The pioneering analysis by Field and Feyn- 
man [U[7] provides a simple and lucid picture to under- 
stand the fragmentation for a quark into mesons based on 
a recursive principle, and has stimulated more sophisti- 
cated models like the string fragmentation model, imple- 
mented in Monte Carlo generation programs yQ , as in the 
Lund model (B|. The Field-Feynman model is successful 
to parameterize the fragmentation functions of the pion 
and kaon with only two parameters. However, no distinc- 
tion between D^ (z) and (z) was made in the origi- 
nal model [B] , for the sake to reduce independent parame- 
ters. With the recent progress of clear distinction between 
D% + (z) and D^ + (z), it becomes mature and necessary to 
consider the flavor structure of the unfavored fragmen- 
tation functions of mesons in the Field-Feynman model. 
The purpose of this paper is, for the first time, to make 
a distinction between (z) and D^ + (z) in the Field- 
Feynman model. It will be shown that this can be realized 
by adopting the flavor dependent parameters for the exci- 
tation of sea quark-antiquark pairs by the initial quarks. 
We will find that the obtained analytical results for D* , 
.Dj , and D* + with only three parameters, are compati- 
ble with the recent KLC parametrization Eqs. (IIIMII based 
on experimental results. We will also extend the same 
framework to the fragmentation functions of the kaon, and 
show that we can distinguish between D^ (z), D^ (z), 
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(z), and D*f + (z) without introducing any additional 
parameters. Thus the predictions can be tested by future 
experiments on the fragmentation functions of the kaon. 

The paper is organized as follows. In Sec. 2, we present 
a brief review of the Field-Feynman recursive model for 
the quark to meson fragmentation functions, and make 
a distinction between the parameters that can introduce 
the difference between (z) and (z). In Sec. 3, we 
adjust the parameters based on the KLC parametriza- 
tion and the experimental results. We present the ana- 
lytical results of , , and , and compare them 
with the experimental results. In Sec. 4, we extend the 
framework to the fragmentation functions of the kaon, and 
make predictions of the four different fragmentation func- 
tions Dl <+ [z), D% + ( z ), Df + 0), and D$ + (z). Finally, we 
present a summary in Sec. 5. 



2 A brief review of the recursive model 

The Field-Feynman model [BJ of the meson fragmentation 
is on the basis of a recursive principle, as can be illustrated 
in Fig. ^ The model ansatz is based on the idea that an 
incoming quark "A" combines with an antiquark "B" from 
a quark-antiquark pair U BB" which is produced from the 
color field excited by the incoming quark "A". Then the 
meson "AB" is constructed, and is ranked as 1st primary 
meson. The remaining quark "B" goes on the same way 
as the quark "A" , and combines with "C" from a quark- 
antiquark pair "CC" which is excited by the quark "B" . 
It constructs the primary meson "_BC" , ranked as 2nd 
primary meson. Then the remaining quark "C" goes on 
the same way in the remaining cascade. 



D 



DC 



CB 



BA 




DD 



CC 



BB 



Fig. 1. Illustration of the Field-Feynman recursive model for 
the fragmentation of mesons. The initial quark A combines 
with B to form the rank 1 "primary" meson with configuration 
BA, and the remained quark B goes on the same way as the 
former quark A. 



With f(rj)dr) denoting the probability that the first 
rank 1 meson leaves a fractional momentum ry to the re- 



maining cascade, the function f(rj) is normalized by 



f(z)dz = 1. 



(4) 



Then let F(z)dz be the probability of finding a meson 
(independent of rank) with fractional momentum z within 
dz in a quark jet. So F(z) can be written by the recursive 
principle 



F(z) = f(l-z)+ t^Lf^Fd). 

Jo V V 



(5) 



The first term comes from the ranked 1 primary meson. 
The second term calculates the probability of the produc- 
tion of a higher rank meson recursively. Field and Feyn- 
man 6 found a solution for the above integral equation 
with 

zF(z) = f(l-z), (6) 

by choosing 



f(z) = (d+l)z d . 



(7) 



With Pij being the probability of finding the qjqj pair 
in the quark sea excited by the initial quark i, the nor- 
malization condition imposes that 



(8) 



Here, we add an i index which does not appear in the orig- 
inal model, and we will find that it is useful. The modifica- 
tion of this model is based on the difference between (3 UU 
and (3 SU , which were assumed to be the same in the original 
Field-Feynman model for the sake to reduce independent 
parameters. The original model predicts no distinction be- 
tween and , and the results in |3j show that this 
is not the realistic situation. Adopting SU(2) symmetry 
between u and d quarks 



(3 iu = I3 ld ee f3 h (9) 



we get 



ft s = 1 - 2ft. (10) 

For an initial quark of flavor q, the mean number of frag- 
mented mesons with configuration ab at z is given, in anal- 
ogy of Eq. ©, by 

Pf(z) = 8 aq (3 gb f(l-z)+ [ ^-f(ri)p qc P?{z/rj). (11) 
The mean number of meson states averaged over all quarks 
P <^) = E«?W, (12) 



so that 



P< q >(z) = /V/W(! -*)+ I jfiv)P q cPZ>(z/v)- 

(13) 
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Comparing with Eq. one yields 



3 The fragmentation functions of the pion 



Pf q> {z) = P qa P qb F(z). 
Substituting it back into Eq. I|13|) . one yields 



(14) 



Pf{z) = 6 aq f3 qb f(l -z) + P qa p qb F(z), (15) 



where 



F{z) = F(z) - /(l -z) = a- 1)/(1 - z) 

= (d+l)z- 1 (l-z) d+{ . [W) 

The fragmentation function is 

Dh M=Y. r al P f^ 

ab 

where r h - is the probability of a meson h with configu- 



„0 



ration ab. For example, r*- = 1 and 7T- = ri = 1/2. 

i ta uu act ' 

Combining Eqs. and l|13[>. we obtain 



£>J(z) = ^/(l-2) + S^), 



where 



4 = £ r >- 

6 

a . b 



(18) 
(19) 
(20) 



Then the fragmentation functions of , D^ + , and D™ + 
can be written as 



DZ (z)=(3 u f(l-z) + (f3 u ) 2 F(z), (21) 

D« d + {z) = {p u YF(z), (22) 

cf(2) = (/3 stl ) 2 F(z) = (/3 s ) 2 F(z), (23) 

here we adopt the notation [3 S — f3 su . 

In the original Field-Feynman model, the parameters 
P u and f3 s are identical: (3 U — p s . This means that the 
effect of producing quark-antiquark pairs excited by the 
incident s or u quark is assumed to be the same. This 
leads to the result of (z) = (z). Indeed, the effect 
due to the quark mass in the hadronization process strictly 
speaking cannot be ignored. In other words, the excitation 
of a quark-antiquark pair by incident quarks with different 
masses will be different. This can be taken into account by 
introducing (3 U ^ p s , and thus obtain (z) ^ (z). 
That is, we take into account the difference in the exci- 
tation of quark-antiquark pairs by the initial light-flavor 
(u or d) and strange (s) quarks. Here, 5/7(2) symmetry is 
still assumed, P qu — P q d — (3 q . 



From (|21|l and (|22l) . two relations can be obtained 

/3u(lA-l) 



Df(z)/DZ (z) 



(1 " Pu + Pu/z) ' 



(24) 



DZ {z)-Dl (z)=p u f(l-z). (25) 

The first relation is independent of f(z), so that it is 
useful to set the parameter j3 u . The parametrization re- 
sults of (z) and (z) are used to adjust the 
parameters. We try to set the value of f3 u by plotting 
D J (z)/DZ (z) for the analytical result in this work and 
the parametrization result by KLC, as shown in Fig. [3 
We find that f3 u = 0.46 in the domain of z € (0.3,1). 
We also adopt the option of (3 U = 0.4, which was used in 
the original model , as a comparison. Then plotting the 
function D^ + (z) — (z) in this work and in the KLC 
parametrization, we can get d = 1.45 with j3 u — 0.46, 
as shown in Fig. [21 It is surprising that we can get com- 
patible curves by only adjusting the parameters (3 U and 
d at z > 0.3, as f(z) = (d + l)z d is only an assump- 
tion in the original model. There is no general reason for 
this choice except for the feasibility of solving the integral 
equation J5j) in an exact way. And now, at some degree, 
the KLC parametrization results extracted from experi- 
mental results without any assumption about favored and 
unfavored transitions, support the Field- Feyman assump- 
tion of f(z) as quite right in the domain of z € (0.3, 1). 
Substituting /3 U = 0.46 and d = 1.45 into (ED and (22J, 
we get the analytical results of (z) and DJ + (z), 



DZ (z) = 1.127(0.54 + 0.46/z)(l - z) 



1.45 



Dl (z) = 0.51842z~ 1 (l - z) 



2.45 



(26) 
(27) 

We compare the analytical results of (|26|) and (|27|l with 
the experimental results in Fig. 0] and Fig. [5] It is obvious 
that the analytical results of l|26() and (|27() can well explain 
the experimental results, especially in the domain of z G 
(0.3,1). 

We now look at the unfavored fragmentation func- 
tion D1 (z), which was predicted to be the same as the 
unfavored fragmentation function D^ (z) in the original 
model. We present the experimental results of D™ + (z) 
and the analytic result l|23|l of this work in Fig. H3 with 
the value of (3 S adjusted to be 0.38. That means that 
38% quark-antiquark pairs in the color field excited by 
an s quark are the light-flavor uu or dd pairs. Whereas 
there are 46% uu or dd quark-antiquark pairs excited 
by a light-flavor u or d quark. We should notice that 
(3 SS ^ (3 US , as p„ = 1 - 2(3 SU = 1 - 2/3 s = 0.24 is 
the fraction of ss pairs excited by an s quark, whereas 
p us = 1 — 2p uu — 1 — 2p u = 0.08 is the fraction of ss 
pairs excited by a u or d quark. This is quite reasonable, 
considering that there should be more fraction of ss pairs 
excited by an s quark than that by a light-flavor (u or d) 
quark. Therefore we have 



Df{z) = 0.35378z" 1 (l - z) 



2.45 



(28) 
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x Ratio of Parametrization Result 

0.6 - * 




0.2 0.4 0.6 0.8 1.0 



Fig. 2. The dash curve is the ratio of /D™ + from the 
KLC parametrization 0, and the solid and dot curves are the 
analytical results of this work with two different values of f3 u . 



0.9 




0.0 - 

0.3 0.6 0.9 

z 



Fig. 3. The dash curve is Dl - from the KLC 

parametrization and the solid curve is (3 u (d + 1)(1 — z) d 
(Eq. with f3 n = 0.46 and d = 1.45. 



If we choose (3 S = 0.46, (z) is identical to D% (z), and 
we return back to the original model. 

We present in Fig. [7| our analytical results for a com- 
plete set of fragmentation functions of the pion with dis- 
tinction between , , and , and compare them 
with the KLC parametrization results. We find that the 
two sets of fragmentation functions are compatible with 
each other, especially at z > 0.3. Considering that a large 
uncertainty may exist in the parametrization of the unfa- 
vored fragmentation function (z) at small z, we may 
consider the analytical results in this work as an alter- 
native set of fragmentation functions for possible appli- 
cations, for example, to check the sensitivity on different 
fragmentation functions in the predictions of various pion 
fragmentation processes 0. 




z 



Fig. 4. The circles are the experimental results from [3], and 
the solid curve is the analytical result of this work. 



20 ~ ? 

o Experimental Result of D d (z) 




0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0 
z 

Fig. 5. The circles are the experimental results from [3], and 
the solid curve is the analytical result of this work. 

One notices from Figs. |2] and [3] that the analytical re- 
sults differ by those obtained by KLC in the small z re- 
gion, as also happens for the analytical results and the 
experimental data in Figs. 0] and [SJ This feature is easy 
to be understood, as the sea content of the meson has not 
been taken into account in the Ficld-Fcynman model. The 
small z behaviors of the fragmentation functions are also 
influenced by the sea quark-antiquark pairs inside the me- 
son U|. Thus it is necessary to consider this aspect if one 
wishes to describe the small z behaviors of fragmentation 
functions with a more reasonable picture. 



4 The fragmentation functions of the kaon 

We now apply the same framework to the fragmentation 
functions of the kaon. Following the same procedures, we 
can get 

Df (z) = (3 8 f(l -z) + p.(l- 2p s )F(z), (29) 
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Experimental Result of D s (z) 
Analytical Result with p s =0.46 
=0.42 

Analytical Result with p s =0.38 



ry-rr—n- 



0.2 



0.4 



0.6 



Df{z) = 0.09016z- 1 (l-z) 2 45 . (36) 

These analytical results of the fragmentation functions 
of the kaon are plotted in Fig. [S] with four independent 
curves: Df + (z), Df + {z), Df + (z), and Df + (z). There 
should be only three independent curves without the dif- 
ference between (3 U and (3 S , e.g., Df (z) equals to Df (z) 
as in the original model. In principle, the flavor structure 
of the kaon fragmentation functions is more complicated 
than that of the pion, as revealed by the Field-Feynman 
model. However, some available parametrizations [1011111 
I12U13| usually make only distinction between the favored 
fragmentation functions, which are related to the valence 
quarks of the kaon, and the unfavored fragmentation func- 
tions, which are related to the light-flavor sea quarks of 
the kaon. 




Fig. 7. The fragmentation functions of the pion with distinc- 
tion between DZ , D% , and , The thick curves are the 
analytical results of this work, and the thin curves are the KLC 
parametrization results 

Df\z) = (1- 2A,)/(1 - z) +/3„(1 - 2(3 u )F(z), (30) 

Df + (z)=f3 s (l-2f3 s )F(z), (31) 

Df + (z) = Df(z) = Df(z) = (3 U (1 - 2(3 u )F(z), (32) 

where (3 S = (3 SU , (3 U = (3 UU , and /3 qa = 1 - 2/3, . 

With the parameters (3 U = 0.46, (3 S — 0.38, and d = 
1.45 being the same as for the pion, we get a complete set 
of fragmentation functions of the kaon 

Df + (z) = 0.931(0.76 + 0.24/z)(l-z) 1 45 , (33) 

Df + {z) = 0.196(0.54+ 0.46/z)(l - z) 1 45 , (34) 
Df + (z) = 0.22344z- 1 (l - z) 2A \ (35) 



The effect of introducing (3 S + 1 (3 U might be small in the 
fragmentation functions of the pion, but it is amplified in 
the fragmentation functions of the kaon. It docs not only 
introduce the difference between the unfavored fragmenta- 
tion functions Df + (z) and Df + (z), but also causes mod- 
ification to the favored fragmentation functions Df + (z). 
In Fig. |21 we plot the ratio of the favored fragmentation 
functions: Df + {z)/Df + (z). The curve with (3 S = 0.42 is 
adopted as an alternative option for comparison, and the 
effect of f3 s on the unfavored fragmentation function is il- 
lustrated in Fig.EB The curve with (3 S — 0.46 is the result 
in the original model, where (3 S — (3 U — 0.46. Therefore 
we notice a significant modification to the favored frag- 
mentation function Df (z) in comparison with that in 
the original model. We also plot the prediction of the un- 
favored fragmentation function Df (z), and find that it 
differs from Df (z) significantly. These predictions can 
be tested by future experiments, such as the semi-inclusive 
kaon productions in deep inelastic scattering by the HER- 
MES collaboration and the CLAS collaboration et ai. The 
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z 

Fig. 9. These three curves give the ratio between the favored 
fragmentation functions (z)/D^ (z) in this work. It goes 
back to the original Field-Feynman model when /3 S — Pu = 
0.46. 

parameters used in this work can be also constrained by 
further detailed studies, or the framework should be ex- 
tended to more general situation with more sophisticated 
aspects taken into account. 
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Fig. 10. These three curves show the unfavored fragmentation 
function D# (z) with different (3 S in this work. It goes back 
to the original Field-Feynman model when f3 s = (5 U = 0.46. 



strange (s) quarks, we can distinguish between , , 

and D™ + in the Field-Feynman model with only three pa- 
rameters. The analytical results obtained are compatible 
with available empirical results, and can be served as an 
alternative set of completed fragmentation functions of 
the pion. We also extended the same framework to the 
kaon, and provided a new set of fragmentation functions of 
the kaon with clear distinction between (z), D^f + (z), 

Djf + (z), and Dff + (z). Thus this work provides the frag- 
mentation functions of the pion and kaon with a same set 
of only three parameters. The predictions can be tested 
by future experiments on the fragmentation functions of 
the pion and kaon. 
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5 Summary 

In summary, we analyzed the fragmentation functions of 
both the pion and kaon in the Field-Feynman recursive 
model with significant modification. By taking into ac- 
count the difference between the excitation of sea quark- 
antiquark pairs by the initial light-flavor (u or d) and 



